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Abstract. Group-theoretical approach is applied to study behavior of lossless two- 
level atoms in a standing-wave laser field. Due to the recoil effect, the internal and 
external atomic degrees of freedom become coupled. The internal dynamics is described 
quantum mechanically in terms of the SU(2) group parameters. The evolution operator 
is found in an explicit way after solving a single ODE for one of the group parameters. 
The translational motion in a standing wave is governed by the classical Hamilton 
equations which are coupled to the SU(2) group equations. It is shown that the full 
set of equations may be chaotic in some ranges of the control parameters and initial 
conditions. It means physically that there are regimes of motion with chaotic center-of- 
mass motion and irregular internal dynamics. It is established that the chaotic regime 
is specified by the character of oscillations of the group parameter characterizing the 
mean interaction energy between the atom and the laser field. It is shown that the 
effect of chaotic walking can be observed in a real experiment with cold atoms crossing 
a standing-wave laser field. 
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1. Introduction 

In quantum physics the unitary time evolution of a driven quantum system is described 
by the evolution operator equation 



where Uh.{t,to) is a time evolution operator, H is a Hamiltonian and h(t) is a vector- 
function of the system's control parameters. From the abstract point of view, the 
evolution equation ([T]) can be regarded as a differential equation on the group of 
dynamical symmetry. By dynamical symmetry we mean simply that the Hamiltonian 
can be expressed as a linear combination of operators belonging to a finite-dimensional 
Lie algebra with n basic elements. The parameters, gj. , k = 1, 2, . . . , n, of the respective 
Lie group satisfy to a set of n first-order ordinary differential nonlinear equations which 
depend only on the structure of the algebra and on c-number coefficients of the system's 
Hamiltonian or the other governing operator [U El El II]- Thus, the dynamical group 
itself may be considered dynamical system. 

The dynamical-symmetry and Lie-algebraic approach has been successfully applied 
to describe the time evolution of numerous physical systems in different disciplines 
extending from classical mechanics [5], classical optics [HI CO El |9] and quantum mechanics 
[T0| [TT] IT2| [T3l [H] to physics of neutrino oscillations [15J . As to study of dynamics of 
laser driven atoms, this approach has been applied to get Lie algebraic solution of the 
Bloch equations [161 E] • 

The evolution of an isolated quantum system is regular, and the overlap of any 
two different quantum state vectors is a constant in course of time. All the expectation 
values of the quantum variables evolve in a quasiperiodic way at most. It does no matter 
how complicated is a dynamical symmetry of the quantum system under consideration 
and the corresponding Lie algebra. On the other hand, it is well known that even 
simple classical systems may be unstable and demonstrate chaotic behavior [191 120] - 
Classical instability is usually defined as an exponential separation of two initially close 
trajectories in time with an asymptotic rate given by the maximal Lyapunov exponent 
A. Such a behavior is possible because of the continuity of the classical phase space 
where the system's states can be arbitrary close to each other. The trajectory concept 
is absent in quantum mechanics, and the quantum phase space is not continuous due to 
the Heisenberg uncertainty principle. Perfectly isolated quantum systems are unitary, 
and there can be no chaos in the sense of exponential instability even if their classical 
limits are chaotic. What is usually understood under quantum chaos is, in essence, the 
special features of the quantum unitary evolution of the system under consideration (no 
matter how complicated the evolution is) in the region of its control parameter values 
and/or initial conditions at which its classical analogue is chaotic [211 E21 E21 EI] . In fact, 
it is not a special quantum problem. Any type of propagating waves (electromagnetic, 
sound or others), satisfying to a linear wave equation (which is an analogue of the 
Scrodinger equation), has the same property. Wave chaos is the special features of the 
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wave field in the region of control parameters and/or initial conditions at which its ray 
analogue is chaotic [25| [26] . Thus, the quantum (wave) chaos problem is partly the 
problem of quantum (wave)-classical (ray) correspondence. 

Let us describe briefly the interconnection between the dynamical symmetries and 
dynamical chaos in physics of the atom-field interaction. 

(i) The simplest problem is dynamics of a two-level atom at rest in an external laser 
field. From the dynamical symmetry point of view, the 577(2) group, generated by 
the corresponding Hamiltonian, is driven by an external force that is not considered 
to be a dynamical system. It is the case of an external driving. The problem 
has been studied in Ref. [I] where it has been shown that the evolution of atomic 
internal variables in a linearly polarized bichromatic laser field with incommensurate 
frequencies may be very complicated on the Bloch sphere albeit regular. It is 
simply because the dynamics takes place on the two-dimensional surface of the 
Bloch sphere. 

(ii) If we deal with a two-level atom at rest in an ideal cavity and take into account the 
response of the atom to the cavity radiation field, the semiclassical evolution of the 
coupled atom-field system may be chaotic in the sense of exponential sensitivity to 
small variations in initial conditions and/or parameters [271 EHl EHl ESS EH |32j 133] . 
This is the case of so-called dynamical driving [I] when the 577 (2) group, generated 
by the atomic Hamiltonian, is driven by another dynamical system, the field 
one. We deal now with a quantum system, the atom, which is coupled with 
a classical system, the radiation field governed by the Maxwell equations. The 
resulting Maxwell- Scro dinger (Bloch) equations constitute the five-dimensional set 
of nonlinear ordinary differential equations with two integrals of motion, the total 
atom-field energy (which is a constant in the absence of any losses) and the length of 
the Bloch vector. The motion in the phase space takes place on a three-dimensional 
manifold and may be chaotic due to transverse intersection of stable and unstable 
manifolds of hyperbolic points in some ranges of the control parameters, the values 
of the maximal Rabi frequency and the coupling strength [33] . 

(iii) If a two-level atom moves within a standing-wave laser field in an open space, not 
in a cavity, the field may be considered as an external driving but one needs to take 
into account the atomic recoil effect, i.e. changes of the atomic momentum after 
absorption or emission photons. If the atom is not especially cold, we may treat 
its translation degree of freedom classically. It is again the case of the dynamical 
driving with the SU(2) group driven by the dynamical system which is now the 
classical atomic degree of freedom. The governing Hamilton-Scrodinger equations 
constitute the five- dimensional set of nonlinear ordinary differential equations with 
two integrals of motion, the atomic total energy, including the kinetic one, and the 
length of the Bloch vector. The motion in the phase space takes place on a three- 
dimensional manifold and may be chaotic in some ranges of the control parameters, 
the values of the maximal Rabi and atomic recoil frequencies. A number of nonlinear 
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dynamical effects have been analytically and numerically demonstrated with this 
system: chaotic Rabi oscillations [3H E2], Hamiltonian chaotic atomic transport 



and dynamical fractals [361 E3 EE], Levy flights and anomalous diffusion [391 EH]. 
These effects are caused by local instability of the center-of-mass motion in a laser 
field. A set of atomic trajectories under certain conditions becomes exponentially 
sensitive to small variations in initial quantum internal and classical external states 
or/and in the control parameters, mainly, the atom-laser detuning. 

In this paper we consider the physical situation mentioned in the third part of our 
nomenclature to focus at the ultimate reasons of chaotic atomic external and internal 
motion and its connection with the SU(2) dynamical symmetry. 

2. Lie algebraic solution for the evolution operator with the SU(2) 
dynamical symmetry 

In a variety of physical problems SU(2) appears to be a group of dynamical symmetry. 
It is known [H H] that the set of three equations for the SU{2) group parameters 
can be reduced to a single second-order differential equation. The form of this 
governing equation depends on the choice of the basis and its exponential ordering. The 
appropriate choice of parameterization of the dynamical group is especially important 
if we need to solve explicitly the governing equation for a given physical Hamiltonian. 

The Hermitian Hamiltonian of a quantum system with the SU(2) dynamical 
symmetry can be cast in the general form 

H(t) = h (t)R + h*(t)R-+h(t)R + , (2) 

where Rq an d R± are the generators that satisfy the commutation relations 



R— , 



-2R , 



R , R± 



±R ± . (3) 



It is convenient to choose the following noncanonical parameterization of the SU(2) 
group 

t 

V = exp[( 9o - ijh dr)^} exp 9 _it exp s+ fl + . (4) 



Substituting Eq.(@| into Eq.flTJ), one finds the set of differential equations for the 
group parameters that can be reduced to the single equation for the group parameter 
9 = exp(#o/ 2 ) 

d 2 g fdh/dt \dg 2 dg 

-d¥-(— + lh VTt +lhl 9 = 0,g(0) = l,-(0) = 0. (5) 

Once Eq.flH} is solved analytically, all the other parameters in the product (HI) may be 
expressed in terms of the parameter g as follows: 

t t 

9- = 1 AAL^L e xp^-i J h dT^J , = -^exp(i J h dT^J . (6) 

o o 
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9 = 9-1 9- 



(7) 



Then any group element in the parameterization 
complex numbers g and g obeying the condition 



9 I + I 9 



~ |2 



1 . 



can be described by a pair of 



(8) 



It should be noted that all these formulas are valid within any representation and 
within any realization of the SU(2) group. It is well known that the unitary irreducible 
representations of SU(2) are characterized by half-integer and integer numbers j. The 
dimensionality of the jth representation is equal to 2j + 1. In the (2j + 1) -dimensional 
space of representation there is a canonical basis 

\j,m), m = -j, -j + . (9) 

The representation matrix elements in the noncanonical parameterization fll]) are given 
by ® 



rjU) 



exp 



— iwl J hQ{r)di 
o 



(J+Q! 



J 

E 

l=-j 



(j-m'y.(j-m)\ 
(j+m')\(j+m)\ 



1/2 



(10) 



\9 



(9) 



I —ml ^ ~*y — m 



(j-l)\(l-m)\{l-m')\- 

To analyze the dynamics of a two-level quantum system we need the two- 
dimensional representation of the SU(2) group. In this case the generators -R's are 
connected with the familiar Pauli matrices 



Ro 



where 



0"+, cr_ 



1 
-1 



[o-z, o-±\ 



R. 



±26-4 



O- 





1 



1 




,(11) 



(12) 



In this representation the Hamiltonian of a driven two-level system has the form 

Hit) = ^hw a a z + hQ*(t)a^ + hTt{t)a + , 

where Q(t) is a time-dependent function which is, in general, a complex-valued one. The 
temporal evolution of the two-level system is now governed by the equation 



(13) 



d 2 g ^dVt/dt 



dt- v n +lWa 

The evolution matrix in the basis 
1 



^+ | O \*g = , g(0) = 1 , ^(0) = . 



is given by 



|1> 



[>(l/2) 



2 ' 2 



£>>|2> 



-iu a t/2 



J dt 



1 1 

2 '2 





dt 




(14) 
(15) 

(16) 
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3. The SU(2) group— Hamilton equations for a two-level atom moving in a 
standing-wave laser field 

We consider a two-level atom with mass m a and transition frequency u a , moving with 
the momentum P along the axis X in a one-dimensional classical laser standing wave 
with the frequency Uf and the wave vector kf. In the frame, rotating with the frequency 
Uf, the model Hamiltonian is the following: 

P 2 1 

H = + ^ K^a ~ u s )b z - h~Q (<t_ + a + ) cos k f X, (17) 

where Qq is the maximal Rabi frequency which is proportional to the square root of the 
number of photons in the wave. The laser field is assumed to be strong enough, so we 
can treat the field classically. 

In the process of emitting and absorbing photons, atoms not only change their 
internal electronic states but their external translational states change as well due to 
the photon recoil. If the atomic mean momentum is large as compared to the photon 
momentum hkf, one can describe the translational degree of freedom classically. The 
position and momentum of a point-like atom satisfy classical Hamilton equations of 
motion which we represent in the normalized form 

x = u r p, p — — < <7-(t) + o- + (t) > sin x, (18) 

where x = kfX and p = P/hkf are normalized classical atomic center-of-mass position 
and momentum, respectively. The dot denotes differentiation with respect to the 
dimensionless time r = Qyt and u) r = h~k"j / m a Qo <C 1 is the normalized recoil frequency. 
To compute the quantum expectation value < <7_(i) + > we need to use the 

solution for the evolution operator ( JTBI) . Supposing that the atom is initially in the 
ground state |1), we get 

< <r_(t) + a + (t) >= (l\fr(t)U(t)\l) = -(gG* + g*G), (19) 

where we introduce for convenience the new complex-valued variable 

G = — (20) 
cosx 



The internal atomic dynamics is governed by Eq. (1141) that can be rewritten in the 
form of two first-order equations for the complex-valued group parameters g and G. The 
self-consistent set of equations for the coupled external and internal atomic degrees of 
freedom now reads as 

x = u r p, p = (gG*+g*G) sin x,g = iG cosx, G = — iAG+ig cosx, (21) 

where the normalized recoil frequency u r and the atom-field detuning, A = (ojf—uj a )/£l , 



are the control parameters. The six-dimensional dynamical system (jzTj) has two 
independent integrals of motion, the total energy, 

H = —p 2 + (gG* + g*G) cosx - -(GG* - gg*), (22) 
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and the integral, 

I 9 | 2 + I G | 2 = 1 (23) 

reflecting conservation of the norm of the atomic wave function. It is evident from the 
second integral (|23|) that the squared absolute values of the SU(2) group parameters, 

I 2 l 2 

| g | and | G | , have the sense of the probability amplitudes to find the atom in the 
ground and excited states, respectively. 

The equations of motion (|2ip describe the mixed quantum-classical dynamics of a 
two-level atom in a one-dimensional standing-wave laser field. The dynamical SU(2) 
group is responsible for internal atomic dynamics caused by the interaction of the atomic 
electric dipole moment with the strength of the electric component of the field. The 
quantum expectation value of the corresponding interaction energy is given by the 
combination of the SU(2) group parameters ffl9l) . The classical translational degree 
of freedom is described by the Hamilton equations (see the first two equations in the set 
( 121]) ) governed by the interaction energy. In Introduction we called such a situation 
as a dynamical driving when the SU(2) group, generated by the atomic quantum 
Hamiltonian, is driven by another dynamical system, the classical atomic degree of 
freedom. In fact, we deal not with a fully quantum system but with a quantum-classical 
hybrid which is described by the c-number nonlinear dynamical system f l2TT) that may be 
chaotic in the strict sense of this term in some ranges of the control parameters and/or 
initial conditions. 

4. Dynamical chaos in the group-theoretical picture 

Equations ( 12T|) constitute a nonlinear autonomous dynamical system with three degrees 
of freedom and, in general, with the two integrals of motion, ( 122]) and ( |23|) . Thus, the 
dynamical system (12TT) may be chaotic in the sense of exponential sensitivity to small 
variation in initial conditions and/or the control parameters, u r and A. The common 
test to confirm that is to compute the maximum Lyapunov exponent characterizing 
the mean rate of exponential divergence of initially close trajectories which serves as a 
quantitative measure of dynamical chaos jlQl HI] : 

A(Q ,Aq )= lim lln ^f^fK (24) 

t^oo,Aq -+0 t || Aqo|| 

where Aq is the vector in the phase space with the components {Aqj,j = 1, N} and 
the norm ||Aq||. In Eq. (|24"]) . Aqo and Aq(Qo,£) denote the separation between two 
initially adjacent trajectories at the initial moment t — and at time t, respectively, 
Qo is the initial position. If, at least, one of the Lyapunov exponents of the dynamical 
system under question is positive, then trajectories, starting close together in the phase 
space, separate exponentially as time grows. This very sensitive dependence on initial 
conditions is one of the main indicator of dynamical chaos. 

The result of computation of the maximum Lyapunov exponent with the equations 
of motion fl2Tl) at u r = 10~ 3 in dependence on the detuning A and the initial atomic 
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momentum p Q is shown in Fig. [TJ Color in the plot codes the value of the maximum 
Lyapunov exponent A. In white regions in Fig. [T]the values of A are almost zero, and the 
atomic motion is regular in the corresponding ranges of A and po- In shadowed regions 
positive values of A imply unstable motion. The atoms with zero A's either oscillate in 
a regular way in a well of the optical potential or move ballistically over the hills of the 
potential with a regular variation of their velocity. 




0.0 0.2 0.4 0.6 0.8 A 1.0 



Figure 1. Maximum Lyapunov exponent, A, vs the atom-field detuning A (in units 
of the maximal Rabi frequency f2) and the initial atomic momentum po (in units of 
the photon momentum Kkf) at oj r = 10~ 3 . Color codes the values of A. 




Figure 2. Regimes of motion of two-level atoms in a one-dimensional deterministic 
standing-wave laser field. Trajectories in the real space at ui r = 10~ 3 : regular flight 
(RF, A = 0.8, po = 45), chaotic flight (CF, A = 0.2, p = 45), chaotic walking (CW, 
A = 0.2, po = 10) and trapping in a potential well (T, A = —0.2, po = 5). ir is in units 
of the wavelength A / . 

At exact resonance, the equations of motion ( 121~|) become integrable due to an 
additional integral of motion, gG* + g*G = const, and we get A = 0. Thus at A = 0, 
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Figure 3. Plots with 50 atomic trajectories with different values of the initial atomic 
momentum < po < 50 but with the same initial position xq — and the same other 
initial conditions, (a) Real space, (b) Momentum space. 

the center-of-mass motion and the motion in the space of the SU(2) group parameters 
are regular. 




Figure 4. Behavior of the mean atom-field interaction energy g\G\ + 52G2 in the 
regimes of (a) regular oscillations in a well of the optical potential and (b) regular 
flight. 

There are three possible chaotic types of motion of a two-level atom in a one- 
dimensional standing-wave laser field. In dependence on the initial conditions and the 
parameter values atoms may oscillate chaotically in a well of the optical potential, 
move ballistically over the hills of the potential with chaotic variations of their velocity 
or perform a chaotic walking. In the regime of the chaotic walking an atom in a 
deterministic standing-wave field alternates between flying through the standing-wave 
and being trapped in the wells of the optical potential. Moreover, it may change the 
direction of motion in a random-like way [38]. We would like to stress that local 
instability produces chaotic center-of-mass motion in a rigid standing wave without any 
modulation of its parameters in difference from the situation with atoms in a periodically 
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Figure 5. The same as in Fig. @] but in the regimes of (a) chaotic flight and (b) chaotic 
walking. 



kicked optical lattice [121 H31 Hi]- To illustrate different types of motion we plot in Fig. [2] 
four trajectories of the atoms in the real space at u r = 10~ 3 corresponding to a regular 
flight (RF), chaotic flight (CF), chaotic walking (CW) and trapping in a potential well 
(T). 
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Figure 6. Projections of the single trajectory of a trapped atom in the six-dimensional 
phase space on the plane of the complex-valued SU(2) group parameter g — gi + igi 
at (a) r = 100, (b) r = 500 and (c) r = 1000. 
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Figure 7. The same as in Fig. [6] but for a regular flight. 
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Let us estimate the values of the control parameters and the initial conditions under 
which atoms oscillate in the first well of the optical potential, move ballistically or walk 
chaotically. At small detunings, A < 1, the total energy ( 122]) consists of the kinetic 
one, K = u r p 2 /2, and the potential one, U = (gG* + g*G) cosx = (g\Gi + g 2 G 2 ) cosx, 
the sum of which is conserved approximately in course of time. The maximal absolute 
value of the optical potential energy is 1. Let the atom is prepared in the ground state 
|1>, i.e., 9i (t = 0) = 1, g 2 (r = 0) = G^r = 0) = G 2 (r = 0) = and U = 0. If 
Kq > I [/max I = 1> then the atom will move ballistically. This occurs if the initial atomic 
momentum, p , satisfies to the condition p > ^j2/uo r ~ 44 at u r = 10~ 3 . If the initial 
conditions are chosen to give < H = K + Uq < 1, then the atom performs a chaotic 
walking. This occurs at < po < 44. The atom will be trapped in the first well of 
the optical potential if Hq < 0. It is posiible with the initial conditions chosen only if 
A < 0. 

To demonstrate strong dependence of the atomic motion on initial conditions 
we compute 50 trajectories with different values of the initial atomic momentum, 
< Po < 50, but with the same initial position, xq = 0, and the same other initial 
conditions. Figure [3] gives an impressive image of dynamical chaos with atoms in a 
laser field both in the real and momentum spaces. Most of the atoms in this bunch 
(with < po < 44) walks chaotically, changing the direction of motion in course of time. 
Atomic trajectories with close initial conditions diverge in the real one-dimensional space 
in such a way that it is practically impossible to predict their final position after the 
predictability time 

T ^X ln A^)' (25) 

where Ax is the confidence interval and AxrO) is the practically inevitable error in 
measuring the initial atomic position. 

It follows from ( )2T|) that the translational motion is described by the equation for 
a nonlinear physical pendulum with the frequency modulation 

x — 1uj r {g\G\ + g 2 G 2 ) sin a; = 0. (26) 

It is clear that the regime of the center-of-mass motion is specified by the character 
of oscillations of the group parameter, giGi + g 2 G 2 , which has the sense of the mean 
interaction energy between the atom and the laser field (see the integral of motion (|22|) ). 
If the atom is trapped in the first well of the optical potential, its center of mass oscillates 
between the first negative and positive nodes, — tt/2 < x < ir/2. If, in addition, the 
control parameters are chosen in appropriate way, it will oscillate periodically. This case 
is shown in Fig. H]a (A = —0.2, p = 5) with periodic albeit modulated oscillations of 
the quantity giGi + g 2 G 2 . Figure H]b is plotted for another regime of the center-of-mass 
motion, a regular ballistic flight with A = 0.8 and p = 45. The quantity g±Gi + g 2 G 2 
again oscillates periodically but with the modulation period that is equal to the flight 
time between two adjacent nodes of the laser standing wave, Tf ~ rr/urPo ~ 70. 
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Behavior of the group parameter g\G\ + g 2 G 2 is absolutely different in the chaotic 
regimes of motion, CF and CW. In the regime of chaotic ballistic flight (see Fig. [5] a with 
A = 0.2 and po = 45), shallow oscillations of that quantity are interrupted by jumps of 
different amplitudes that occur when the atom crosses each node of the standing wave. 
In the regime of chaotic center-of-mass walking (see Fig. [5]b with A = 0.2 and po = 10), 
oscillations of the quantity g\G\ + g 2 G 2 even more complicated. We may conclude 
that namely the chaotic oscillations of the mean interaction energy between the atom 
and the laser field, g\G\ + g 2 G 2 , in some ranges of the control parameters, u r and A, 
and initial atomic momentum p are responsible for the chaotic center-of-mass motion. 
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Figure 8. The same as in Fig. [6] but for a chaotic flight. 
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Figure 9. The same as in Fig. [6] but for a chaotic walking. 

The equations of motion 02 II) can be recast in the form of the two second-order 
differential equations, the classical one fl26|) . describing the center-of-mass motion, and 
the quantum one 

g + (zA + xtanx)g + g cos 2 x = 0, (27) 

describing the internal atomic dynamic in terms of the complex-valued SU(2) group 
parameter g = gx + ig 2 . In order to illustrate how different may be behavior of the 
quantum degree of freedom of the quantum-classical hybrid, we compute the evolution 
of the real and imaginary parts of g in course of time. The results are shown in Figs. [6]- 
M with different regimes of motion. The plots give projections of the single atomic 
trajectory in the six- dimensional phase space on the plane of the complex-valued SU(2) 
group parameter g at the time moments r = 100, r = 500 and r = 1000. 
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The plots with a periodically oscillating atom in a trap (Fig. [6]) and with a 
regular flight (Fig. [7]) demonstrate the strictly periodic patterns on the g\—gi plane 
with forbidden regions in the center. Internal dynamics of the atoms in the chaotic 
center-of-mass regimes of motion, chaotic flight in Fig. [8] and chaotic walking in Fig. HI 
is much more complicated. In both the cases, the trajectories on the g\-gi visit in course 
of time all the accessible part of the plane with | g \ < 1 . 



5. How to observe chaotic walking of atoms in a real experiment 

In this section we propose the scheme of a real experiment to observe the effect of 
chaotic walking of atoms in a deterministic standing wave described in the previous 
section. A beam of two-level atoms in the z direction crosses a standing-wave laser 
field with optical axis in the x direction (Fig. [TUh). One measures either the atomic 
density on a substrate as in the atom-lithography experiments H7] or the spatial 
atomic distribution as in the atom optics experiments [4*2l |4"3| HI] . In each type of the 
experiments the results are expected to be different in the regimes of chaotic atomic 
walking and regular motion. To switch between the regimes it is enough to vary the 
value of the detuning in the appropriate way. The laser beam has the Gaussian profile 
exp[— (z — zq) 2 /r 2 ] with r being the e~ 2 radius at the laser beam waist. The longitudinal 
velocity of atoms, v z , is much larger than their transversal velocity v x and is supposed 
to be constant. Therefore, the spatial laser profile may be replaced by the temporal one. 
The Hamiltonian ( TT7|) now takes the time- dependent form 

P 2 h 3 
H = ^- + -(u a -u f )6- z -hn exp[-(v z t--r) 2 /r 2 } (<x_ + o+) cos k f X(28) 



with the same dynamical symmetry. Using the same normalization as before, we get 
the equations of motion 

x - tu r n(r) (gG* + g*G) sin x = (29) 

0(r) 



g + 



iA + xt&nx 



0(r) 



g + g[n(r)Y cos 2 x = (30) 



with the time-dependent coefficient CI(t) = exp[— (r — |cr r ) 2 /cr^], where o T = tVLq/v z is 
the normalized characteristic interaction time. 

To be concrete let us take lithium atoms with the relevant transition 2S , 1 / 2 — 2P 3 / 2 , 
the corresponding wavelength A a = 670.7 nm and the recoil frequency z/ rec = 63 KHz. 
With the maximal Rabi frequency f2 /27r ~ 126 MHz and the radius of the laser beam 
r = 0.05 cm one gets u r = 10~ 3 and a T = 400. To simulate a real experiment let us 
consider a beam of 10 4 lithium atoms with the initial Gaussian position and momentum 
distributions (the rms o x = o v = 2, the average values, x = 0, and p = 10) and 
compute their position distribution at a fixed moment of time. In Fig. [10] b we compare 
the atomic position distributions at r = 1000 (z = 200 microns) for the chaotic walking 
at A = 0.2 (bold curve) and the regular motion at A = 1 (dashed curve). The difference 
is evident. In the chaotic regime atoms are distributed more or less homogeneously over 
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Figure 10. (a) Scheme of the proposed experiment on observation of chaotic walking 
(cw) of atoms scattered at a Gaussian standing laser wave, (b) The distributions of 10 4 
lithium atoms at r — 1000 (z = 200 microns) under the conditions of chaotic walking 
at A = 0.2 (bold curve) and regular motion (rm) at A = 1 (dashed curve). 



a large distance of 8 wavelengths along the x-axis whereas in the regime of the regular 
motion they form a few peaks in a much more narrow interval. Thus, we predict that 
under the conditions of chaotic walking there should appear a less contrast and more 
broadened atomic relief as compared to the case of regular motion because a large 
number of atoms are expected to be deposited between the nodes as a result of chaotic 
walking along the standing-wave axis. 

6. Conclusion 

We have studied behavior of lossless two-level atoms in a one-dimensional standing- 
wave laser field in the group-theoretical picture. In this picture we have represented the 
internal quantum atomic dynamics in terms of the dynamical SU (2) group parameters 
and the center-of-mass motion by the classical Hamilton equations. Thus, we have 
modeled the system by a quantum-classical hybrid with coupled quantum and classical 
degrees of freedom. We have derived the corresponding set of the 577(2) group-Hamilton 
equations of motion with, in general, two integrals of motion. This set has been 
numerically shown to be chaotic in some ranges of the control parameters and initial 
conditions. We have found five different regimes of the center-of-mass motion including 
chaotic walking when an atom in an absolutely deterministic standing-wave field may 
change the direction of motion in a random-like way alternating between flying in the 
optical potential and being trapped in its wells. All the regimes have been illustrated 
by the trajectory plots in the real and momentum spaces. It has been established that 
the instability of motion and dynamical chaos are caused by the character of oscillations 
of the group parameter characterizing the mean interaction energy between the atom 
and the laser field. Projections of atomic trajectories in the six-dimensional phase 
space on the plane of the complex-valued SU (2) group parameter g have been shown to 
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form regular and irregular patterns in the regimes of regular and chaotic center-of-mass 
motion, respectively. 

We proposed the scheme of an experiment on the scattering of atomic beams at 
a standing-wave laser field that could directly image chaotic walking of atoms along 
the optical axis. In a real experiment the final spatial distribution can be recorded 
via fluorescence or absorption imaging on a CCD, commonly used methods in atom 
optics experiments yielding information on the number of atoms and the cloud's spatial 
size. The other possibility is a nanofabrication where the atoms after the interaction 
with the standing wave are deposited on a silicon substrate in a high vacuum chamber. 
In this case the spatial distribution can be analyzed with an atomic force microscope. 
The modern tools of atom optics enable to create narrow initial atomic distributions in 
position and momentum, reduce coupling to the environment and technical noise, create 
one-dimensional optical potentials, and to measure spatial and momentum distributions 
with high sensitivity and accuracy [22j H3J HI] . 

The results obtained can be applied to other models of the atom-field interaction as 
well. In particular, relaxation processes in two-level atoms can be described within the 
framework of the 50(3) dynamical-symmetry approach to solving the Bloch equations 
[T7] . Moreover, one may consider by the method developed in this paper the dynamics 
not only of two-level atoms but of three-, four- and multilevel atoms excited by a few 
laser fields at different atomic transitions. If the corresponding model Hamiltonian has 
the SU(2) dynamical symmetry, then one may use the solution obtained in Sec. 2 that 
is valid for any representation of the SU (2) group. 

The model considered can be generalized to the case with two- level atoms inside a 
high-quality cavity with a quantized field. In the rotating wave approximation the state 
space of the corresponding Jaynes-Cummings model splits up into an infinite class of two- 
dimensional non-communicating subspaces each of which being labeled by eigenvalues 
of the Casimir operator. The system evolves in such a way that transitions between the 
subspaces with different eigenvalues are forbidden. The solution for the time-evolution 
operator in each of these subspaces is given by the matrix (1151) with the group parameter 
satisfying to the equation similar to (JSJ). The resulting equations of motion for the 
coupled atom-field system are expected to constitute an infinite-dimensional set of the 
type ( I2TT) with the group equation ( flBl) acting in each of the subspaces labeled by its 
own eigenvalue. This set is expected to admit very different regimes of motion including 
chaotic ones. 
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